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Abstract 

This paper proposes a robust control method based on sliding mode design for two-level quantum 
systems with bounded uncertainties. An eigenstate of the two-level quantum system is identified as 
a sliding mode. The objective is to design a control law to steer the system's state into the sliding 
mode domain and then maintain it in that domain when bounded uncertainties exist in the system 
Hamiltonian. We propose a controller design method using the Lyapunov methodology and periodic 
projective measurements. In particular, we give conditions for designing such a control law, which can 
guarantee the desired robustness in the presence of the uncertainties. The sliding mode control method 
has potential applications to quantum information processing with uncertainties. 

Index Terms 

quantum control, sliding mode control, bounded uncertainty, periodic projective measurement, 
Lyapunov methodology. 

I. INTRODUCTION 

The manipulation and control of quantum systems is becoming an important task in many 
fields JD-El, such as atomic physics flU, molecular chemistry [[51 and quantum information J6]|. 
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It is desirable to develop quantum control theory in a systematic way in order to adapt it to the 
development of quantum technology 0. Several useful tools from classical control theory have 
been introduced to the control analysis and design of quantum systems. For example, optimal 
control theory has been used to assist in control design for closed and dissipative quantum systems 
lf8l- lTT4ll . A learning control method has been presented for guiding the control of chemical 
reactions 0, lfT5l . Quantum feedback control approaches including measurement-based feedback 
and coherent feedback have been used to improve performance for several classes of tasks such 
as preparing quantum states, quantum error correction, controlling quantum entanglement |[T6ll - 
11301 . Robust control tools have been introduced to enhance the robustness of quantum feedback 
networks and linear quantum stochastic systems OTTl . Il32ll . 

Although some progress has been made, more research effort is necessary in controlling quan- 
tum phenomena. In particular, the robustness of quantum control systems has been recognized as 
a key issue in developing practical quantum technology If3~3~l - ll3~5l . In this paper, we focus on the 
robustness problem for quantum control systems. In ||32|| . James and co-workers have formulated 
and solved a quantum robust control problem using the H°° method for linear quantum stochastic 
systems. Here, we develop a variable structure control approach with sliding modes to enhance 
the robustness of quantum systems. The variable structure control strategy is a widely used 
design method in classical control theory and industrial applications where one can change the 
controller structure according to a specified switching logic in order to obtain desired closed-loop 
properties 11361 . Il3~7i In |[38l and ll3~9ll , Dong and Petersen have formulated and solved a variable 
structure control problem for the control of quantum systems. However, the results in |[38l only 
involve open-loop control design using an idea of changing controller structures and do not 
consider the robustness which can be obtained through sliding mode control. Ref. Il38l and Ref. 
[40J have briefly discussed the possible application of sliding mode control to quantum systems. 
In ipTfll. two approaches based on sliding mode design have been proposed for the control of 
quantum systems and potential applications of sliding mode control to quantum information 
processing have been presented. Following these results, this paper formally presents a sliding 
mode control method for two-level quantum systems to deal with bounded uncertainties in the 
system Hamiltonian [421 - In particular, we propose two approaches of designing the measurement 
period for different situations which are dependent on the bound on the uncertainties and the 
allowed probability of failure. 



Variable structure control design with sliding modes generally includes two main steps: se- 
lecting a sliding surface (sliding mode) and controlling the system to and maintaining it in this 
sliding surface. Being in the sliding surface guarantees that the quantum system has the desired 
dynamics. We will select an eigenstate of the free Hamiltonian of the controlled quantum system 
as a sliding mode. In the second step, direct feedback control is not directly applicable since 
we generally cannot acquire state information without destroying the quantum system's state. 
Hence, we propose a new method to accomplish this task, which is based on the Lyapunov 
methodology and periodic projective measurements. The Lyapunov methodology is a powerful 
tool for designing control laws in classical control theory and has also been applied to quantum 
control problems [|43l - ll48l . Most existing results on Lyapunov control of quantum systems focus 
on designing a control law to ensure that the controlled quantum system's state asymptotically 
converges to the target state. The existing Lyapunov design methods in quantum control rely 
on perfect knowledge of the initial quantum states and system Hamiltonian. In our approach, 
once the Lyapunov control steers the quantum system into a sliding mode domain, we make 
a projective measurement on the system. Hence, the Lyapunov design method can tolerate 
small drifts (uncertainties) when carrying out our control tasks, which will be demonstrated 
by simulation in Section II.C. Periodic projective measurements are employed to maintain the 
system's state in the sliding mode domain when uncertainties exist in the system Hamiltonian. 
If the measurement period is small enough and the initial state is an eigenstate, the frequent 
measurements make the system collapse back to the initial state. This is related to the quantum 
Zeno effect (for details, see 11491 , ll50l and j5T]|). In contrast to the quantum Zeno effect, our 
objective is to design a measurement period which is as large as possible. The framework of the 
proposed method involves unitary control (Lyapunov control) and projective measurement. In 
this sense, it is similar to the discrete-time quantum feedback stabilization problem in ll52ll and 
ll53l . However, these papers do not consider possible uncertainties in the system Hamiltonian and 
use generalized measurements rather than periodic projective measurements. The main feature 
of the proposed method is that the control law can guarantee control performance when bounded 
uncertainties exist in the system Hamiltonian. 

This paper is organized as follows. Section |n] introduces a quantum control model, defines the 
sliding mode and formulates the control problem. In Section Unl we present a sliding mode control 
method based on the Lyapunov methodology and periodic projective measurements for two-level 



quantum systems with bounded uncertainties. Using the known information about uncertainties 
(e.g., the uncertainty bound and type of uncertainties), we propose two approaches (i.e., Eqs. 
<fT3T ) and (fl4l) ) for designing the measurement period to guarantee the control performance. An 
illustrative example is presented to demonstrate the proposed method. The detailed proofs of the 
main theorems are presented in Section [IV] Conclusions are given in Section |V] 

II. Sliding modes and problem formulation 

In this section, we first introduce a two-level quantum control model. Then a sliding mode is 
defined using an eigenstate. Finally the control problem considered in this paper is formulated. 

A. Quantum Control Model 

In this paper, we focus on two-level pure-state quantum systems. The quantum state can be 
represented by a two-dimensional unit vector \\ff) in a Hilbert space Jf?. Since the global phase 
of a quantum state has no observable physical effect, we do not consider the effect of global 
phase. If we denote the Pauli matrices o = {a x , o y , a z ) as follows: 

we may select the free Hamiltonian of the two-level quantum system as Hq = I z = jO z . Its two 
eigenstates are denoted as |0) and |1). To control a quantum system, we introduce the following 
control Hamiltonian H u = Y,k u k(t)Hk, where w&(?) G R and {H{\ is a set of time-independent 
Hamiltonians. For simplicity, the control Hamiltonian for two-level systems can be written as 

H u = u x (t)I x + u y (t)Iy + u z (t)I z , where 

i i fo i\ i i (o -i\ 

The controlled dynamical equation can be described as (we have assumed h = 1 by using 
atomic units in this paper) 

i\y(t)) =#o|y(0> +Lk=x, y ,z^(t)h\v(t)), 
\ ¥ (t = 0)} = \Y }. 



This control problem is converted into the following problem: given an initial state and a target 
state, find a set of controls {uk(t)} in © to drive the controlled system from the initial state to 
the target state. 

In practical applications, we often use the density operator (or density matrix) p to describe 
the quantum state of a quantum system. For a pure state | iff) , the corresponding density operator 
is p = \ v|- For a two-level quantum system, the state p can be represented in terms of the 
Bloch vector r = (x,y,z) = (tr{po x },tv{po y },tr{po z }): 

p = I(/+r-<7). (4) 

The evolution equation of p can be written as 

p = -i[H,p] (5) 

where [A, B] =AB — BA and H is the total system Hamiltonian. 

After we represent the state p with the Bloch vector, the pure states of a two-level quantum sys- 
tem correspond to the surface of the Bloch sphere, where (x,y,z) = (sin 6 cos (p, sin 6 sin <p, cos 6), 
G [0,4 <p G [0,2tt]. An arbitrary pure state |y/) for a two-level quantum system can be 
represented as 

| V A)=cos||0) + ^sin||l). (6) 

B. Sliding Modes 

Sliding modes play an important role in variable structure control [36]. Usually, the sliding 
mode is constructed so that the system has desired dynamics in the sliding surface. For a quantum 
control problem, a sliding mode may be represented as a functional of the state \\ff) and the 
Hamiltonian H; i.e., = 0. For example, an eigenstate of the free Hamiltonian Hq 

(i.e., Ho\(j)j) =Xj\tyj) where Xj is one eigenvalue of Hq) can be selected as a sliding mode. We 
can define S(\\j/),H) = 1 — \ (\j/\(f)j}\ 2 = 0. If the initial state |y/b) is in the sliding mode; i.e., 
S(|ty/o),#) = 1 — | ( V/b 1 0y ) | 2 = 0, we can easily prove that the quantum system will maintain its 
state in this surface under only the action of the free Hamiltonian Hq. In fact, | = e~ !// ° ? |i//b), 
and we have 

S(\ ¥ (t),H) = 1 - imMM 2 = 1 - \(Wo\e iH ^j)\ 2 

= 1 - |(v/ |f-)^f = 1 - \(Wo\(l>j)\ 2 \e ajt \ 2 ( 7 ) 
= 0. 



That is, an eigenstate of Ho can be identified as a sliding mode. For two-level quantum systems, 
we may select either |0) or |1) as a sliding mode. Without loss of generality, we identify the 
eigenstate |0) of a two-level quantum system as the sliding mode in this paper. 

C. Problem Formulation 

In Section II.B, we have identified an eigenstate |0) as a sliding mode. This means that if a 
quantum system is driven into the sliding mode, its state will be maintained in the sliding surface 
under the action of the free Hamiltonian. However, in practical applications, it is inevitable that 
there exist noises and uncertainties. In this paper, we suppose that the uncertainties can be 
approximately described as perturbations in the Hamiltonian. That is, the uncertainties can be 
denoted as = £ x (t)I x + £ y (t)I y + £ z (t)I z . The unitary errors in |33] belong to this class of 
uncertainties and uncertainties in one-qubit (one quantum bit) gate also correspond to this class 
of uncertainties PTI . For a spin system in solid-state nuclear magnetic resonance (NMR), external 
noisy magnetic fields and unwanted coupling with other spins may lead to uncertainties in this 
class. Further, we assume the uncertainties are bounded; i.e., 



When £ = 0, = 0. That is, there exist no uncertainties, which is trivial for our problem. 
Hence, in the following we assume e > 0. An important advantage of classical sliding mode 
control is its robustness. Our main motivation for introducing sliding mode control to quantum 
systems is to deal with uncertainties. We further suppose that the corresponding system without 
uncertainties is completely controllable and arbitrary unitary control operations can be generated. 
This assumption can be guaranteed for a two-level quantum system if we can realize arbitrary 
rotations along the z-axis and £-axis (C, =x ox y) (e.g., see ll54l for details). 

The control problem under consideration is stated as follows: design a control law to drive 
and then maintain the quantum system's state in a sliding mode domain even when bounded 
uncertainties exist in the system Hamiltonian. Here a sliding mode domain may be defined as 
® = { I W) '■ I (0| y) | 2 > 1 — Po, < po < 1 }, where po is a given constant. Here we assume po ^ 0, 1 
since the case po = only occurs in the sliding mode surface and the case po = 1 is always 
true. Hence, the two cases with po = and po = 1 are trivial for our problem. The definition of 
the sliding mode domain implies that the system has a probability of at most po (which we call 




(8) 



the probability of failure) to collapse out of S> when making a measurement. This behavior is 
quite different from that which occurs in traditional sliding mode control. Hence, we expect that 
our control laws will guarantee that the system's state remains in S except that a measurement 
operation may take it away from S) with a small probability (not greater than po). The control 
problem considered in this paper includes three main subtasks: (i) for any initial state (assumed 
to be known), design a control law to drive the system's state into a defined sliding mode domain 
S; (ii) design a control law to maintain the system's state in S; (iii) design a control law to drive 
the system's state back to S> if a measurement operation takes it away from S. For simplicity, 
we suppose that there exist no uncertainties during the control processes (i) and (iii). 

III. Sliding mode control based on Lyapunov methods and projective 

MEASUREMENTS 

A. General Method 

The first task is to design a control law to drive the controlled system to the chosen sliding 
mode domain S. Lyapunov-based methods are widely used to accomplish this task in traditional 
sliding mode control. If the gradient of a Lyapunov function is negative in the neighborhood of 
the sliding surface, then the controlled system's state will be attracted to and maintained in Si. 
The Lyapunov methodology has also been used to design control laws for quantum systems [|43l - 
[48]. However, these existing results do not consider the issue of robustness against uncertainties. 
Since the measurement of a quantum system will inevitably destroy the measured state, most 
existing results on Lyapunov-based control for quantum systems in fact use a feedback design 
to construct an open-loop control. That is, Lyapunov-based control can be used to first design 
a feedback law which is then used to find the open-loop control by simulating the closed-loop 
system. Then the control can be applied to the quantum system in an open-loop way. Hence, 
the traditional sliding mode control methods using Lyapunov control cannot be directly applied 
to our problem. 

Although quantum measurement often has deleterious effects in quantum control tasks, recent 
results have shown that it can be combined with unitary transformations to complete some 
quantum manipulation tasks and enhance the capability of quantum control HOl . Il55l - ll59l . For 
example, Vilela Mendes and Man'ko iffOll showed nonunitarily controllable systems can be made 
controllable by using "measurement plus evolution". Quantum measurement can be used as a 



control tool as well as a method of information acquisition. It is worth mentioning that the effect 
of measurement on a quantum system as a control tool can be achieved through the interaction 
between the system and measurement apparatus. In this paper, we will combine the Lyapunov 
methodology and projective measurements (with the measurement operator o z ) to accomplish 
the sliding mode control task for two-level quantum systems. The projective measurement with 
o z on a two-level system makes the system's state collapse into |0) (corresponding to eigenvalue 
1 of o z ) or 1 1) (corresponding to eigenvalue —1 of o z ). 
The steps in the control algorithm are as follows (see Fig. 1): 

1) Select an eigenstate |0) of Hq as a sliding mode S(\\j/),H) = 0, and define the sliding mode 
domain as @= : |(0|v/)| 2 > 1 - po}- 

2) For a known initial state |y/b), construct a Lyapunov function V(|y/b),5) to find the control 
law that can drive |y/o) into the sliding mode S. 

3) For a specified probability of failure po and V(\\j/o),S), construct the control period 7b so 
that the control law can drive the system's state into @ in a time period Tq. 

4) For an initial condition which is another eigenstate |1), design a Lyapunov function 
V(|l),5) and construct the period T\ by using a similar method to that in 3). 

5) According to po and £, design the period T for periodic projective measurements. 

6) Use the designed control law to drive the system's state into @ in Tq and make a projective 
measurement at t = Tq. Then repeat the following operations: make periodic projective 
measurements with the period T to maintain the system's state in *3>\ if the measurement 
result corresponds to |1), we use the corresponding control law to drive the state back into 
0. 

From the above control algorithm, we see that the design of Lyapunov functions and the 
selection of the period T for the projective measurements are the two most important tasks. To 
design a control law for quantum systems, several Lyapunov functions have been constructed, 
such as state distance-based and average value-based approaches ll43l - H8l . Here we select a 
function based on the Hilbert-Schmidt distance between a state |y/) and the sliding mode state 
as a Lyapunov function [|46ll . 11481 : i.e., 

V(\ ¥ ),S)= l -(l-\(<j) J \ V }\ 2 ). 



■ Lyapiniov ^ Me asurcment 



uncertainties 




time t 



Fig. 1. The sliding mode control scheme for a two-level quantum system based on Lyapunov methods and periodic projective 
measurements. In this figure, "Lyapunov", "Measurement" and "uncertainties" mean the evolution process of the quantum system 
under the Lyapunov control law, the projective measurement and uncertainties in the system Hamiltonian, respectively. 



It is clear that V > 0. The first-order time derivative of V is 

V =-Zk=x,y,z^mh)^j\h\V)} 

= - E* ,,, u k \ < w\ fc> 13 [e iZ ^ {Hh\ ¥)], 
where 3[a + bi] = b (a,b E R) and Zc denotes the argument of a complex number c. To ensure 
V < 0, we choose the control laws as in ||48l : 

u k = K k f k {Z[e lA ^H<Pj\I k \w)]), (k = x,y,z) (10) 

where K k > may be used to adjust the control amplitude and /(•) satisfies xf(x) > 0. Define 
Z{y\4>j) = 0° when (yr\ty) = 0. 

When one employs a Lyapunov methodology to design a control law, LaSalle's invariance 
principle is a useful tool to analyze its convergence. That is, if x(t) = g(x(t)) is an autonomous 
dynamical system with phase space Q. and V(x) is a Lyapunov function on Q. satisfying V(x) > 
for all x ^ xq and V(x(t)) < 0, any bounded solution converges to the invariant set E = 
{x\V(x(t)) = 0} as ?-)• +°° (for details, see ll60l ). For two-level quantum systems, LaSalle's 
invariance principle can guarantee that the quantum state converges to the sliding mode |0) 
under the control law in (flOl) (for details, see Il48l0 . The convergence is asymptotic. Hence, we 



make a projective measurement with the measurement operator o z when we apply the Lyapunov 
control to the system for Tq (corresponding to the initial condition |Vo}) or T\ (corresponding 
to the initial condition 1 1)), which will drive the system into |0) with a probability not less than 
1 -Po- 

Another important task is to design the measurement period T. We can estimate a bound 
according to the bound e on the uncertainties and the allowed probability of failure p$. Then, 
we construct a period T to guarantee control performance according to the estimated bound. An 
extreme case is T — > 0. That is, after the quantum system's state is driven into the sliding mode, 
we make frequent measurements. This corresponds to the quantum Zeno effect ll50ll . which is the 
inhibition of transitions between quantum states by frequent measurement of the state (see, e.g., 
||49l and 11501 ). Frequent measurements (i.e., 7—^0) can guarantee that the state is maintained in 
the sliding mode in spite of the presence of uncertainties. However, it is usually a difficult task 
to make such frequent measurements. We may conclude that the smaller T is, the bigger the 
cost of accomplishing the periodic measurements becomes. Hence, in contrast to the quantum 
Zeno effect, we wish to design a measurement period T as large as possible. In the following 
subsection, we will propose two approaches of designing T for different situations. 

B. The Design of the Measurement Period T 

We select the sliding mode as S(\\ff),H) = 1 — |(v|0)| = 0- ^ there exist no uncertainties and 
we have driven the system's state to the sliding mode at time to, it will be maintained in the 
sliding mode using only the free Hamiltonian Hq; i.e., S(|i/A( f>fo )),i/o) = 0. That is, if the quantum 
system's state is driven into the sliding mode, it will evolve in the sliding surface. However, in 
practical applications, some uncertainties are unavoidable, which may drive the system's state 
away from the sliding mode. We wish to design a control law to ensure the desired robustness 
in the presence of uncertainties. Assume that the state at time t is p f . If we make measurements 
on this system, the probability p that it will collapse into |1) (the probability of failure) is 

p = <l| Pf |l) = iz* J (11) 

where z t = \x(p t o z ). We have assumed that the possible uncertainties can be described by — 
£x(t)I x + £y(t)I y +£ z (t)I z , where unknown e x (t), e y (t) and e z (t) satisfy ^ e}{t) + e^it) + e}{t) < e. 
We now give detailed discussions to design the measurement period T for possible uncertainties. 



First we consider a special case = e(t)I z (\£(t)\ < £). This case corresponds to phase-flip 
type bounded uncertainties. For any H\ = £,/, (where |e z | < £), if S(\\yo),H) = 0, we have 

S(\y(t)),H) = l-\( ¥ (t)\0}\ 2 

= l-\{\j/ \e'^+ £ ^ t \0)\ 2 



l-\(¥o\0)\ 2 \e 



(12) 



= 0. 

This type of uncertainty does not drive the system's state away from the sliding mode. Hence 
we ignore this type of uncertainty in our method. 



Now we consider the unknown uncertainties // A = £ x (t)I x + £ y (t)I y (where J £ 2 {t) +£y(t) < £) 
and have the following theorem. 

Theorem 1: For a two-level quantum system with the initial state | yy(0) ) = |0) at the time 
t = 0, the system evolves to \y{t)) under the action of H(t) = L + £ x {t)I x + £ y (t)I y (where 



£ 2 (t)+£ 2 (t) < £ and £ > 0). If t e [0,T(% where 

(i) arccos(l — 2po) 



£ 

|2 



(13) 



the system's state will remain in @ = : |(0|i//}| > 1 — po} (where < po < 1). When one 
makes a projective measurement with the measurement operator o z at the time t, the probability 
of failure p = \ (l\\j/(t)}\ 2 is not greater than p$. 

Using Theorem 1, we may try to maintain the system's state in Qi (i.e., the subtask (ii)) by 
implementing periodic projective measurements with the measurement period T = If we 
have more knowledge about the uncertainties, it is possible to improve the measurement period 
Now assume that the uncertainty is = £(t)I^ (£ =x or y) and po G (0, f^j]- We have 
the following theorem. 

Theorem 2: For a two-level quantum system with the initial state |y(0)) = |0) at the time 
t = 0, the system evolves to under the action of H(t) = I z -\-£(t)Ir (where £ = x or y, 

\£(t) \ < £ and £ > 0). If po G (0,^] and t G [0,7^], where 



, (2 ) _ arccos[l-2(l + ^)p ] 



(14) 



VT+e 2 

the system's state will remain in 3 = {\\ff) : |(0|i//)| 2 > 1 — po} (where < po < 1). When one 
makes a projective measurement with the measurement operator o z at the time t, the probability 
of failure p = \ {l\y{t))\ 2 is not greater than pq. 



Remark 1: The proofs of Theorem 1 and Theorem 2 will be presented in Section |IVj 
The two theorems mean the following fact. For a two-level quantum system with unknown 
uncertainties H A = £ x (t)I x + £ y (t)I y (where J~^(t) + £ y (t) < £), if its initial state is in the sliding 
mode |0), we can ensure that the probability of failure is not greater than a given constant 
Po (0 < po < 1) through implementing periodic projective measurements with the measurement 
period T = T^ 1 ) using (TT3T) . Further, if we know that po and £ satisfy the relationship < 

2 

Po < jrgi an d there exists only one type of uncertainty (i.e., Ha = e(t)I x or Ha = £(t)I y , where 
l e (0l — £ )' we can design a measurement period T = T^> using (fl4l) which is larger than 
T^ 1 ). The proof of Theorem 2 also shows that T^> is an optimal measurement period. This 
measurement period will guarantee the required robustness. It is easy to prove the relationship 
j( 2 ) > fW f or arbitrary po e (0, j^gz]- The detailed proof will be presented in the Appendix. 
Based on the above analysis, the selection rule for T is summarized in Table I. Moreover, from 
(fT3l) and (fl4l> . it is clear that for a constant £, r^ 1 ) -)> and ->■ when p ->■ 0. That is, 
for a given bound £ on the uncertainties, if we expect to guarantee the probability of failure 
po — > 0, it requires us to implement frequent measurements such that the measurement period 
T — > 0. Another special case is £ — > +oo 5 which leads to — > and — > 0. That is, to deal 
with very large uncertainties, we need to make frequent measurements (T — > 0) to guarantee the 
desired robustness. From (TT3T) . we also know that for a given po, monotonically decreases 
with increasing £. This means that we need to employ a smaller measurement period to deal 
with uncertainties with a larger bound £. 



Type of uncertainties 


H A = £ x (t)l x +£y{t)ly 


H A =£(t)l c (£=*or y) 


Allowed probability of failure po 


0<p <\ 


0<Po< ^ 


iw <P0<1 


The measurement period T 


T = 


T = T^> 





TABLE I 

A SUMMARY ON THE SELECTION RULE OF THE MEASUREMENT PERIOD T 



C. An Illustrative Example 

Now we present an illustrative example to demonstrate the proposed method. Assume po = 
0.01. Consider two cases: (a) £ = 0.02; (b) £ = 0.2. For simplicity, we assume |ty/o) = 



D.01 0.02 0.03 0.04 0.05 0.06 0.07 

t 



Fig. 2. The probability of |0) under Lyapunov control. 

Hence, Tq = T\. We first design the control and T\ using (flOl) . Here, we consider control only 
using H u = \u(t)o y . Using fll), we select u{t) = K(3 [e ! ' z ^MI°> (0|o>,| y(*))D and K = 100. Let 
the time stepsize be given by 8t — 10~ 4 . We can obtain the probability curve of |0) shown 
in Fig. [21 the control value shown in Fig. [3] and T\ = 0.060. For £ = 0.02, we can design the 
measurement period T = = 10.017 using (fl"3l) . For £ = 0.2, we can design the measurement 
period T = = 1.002 using CL1- Since p' = = 3-8 x 10" 2 > p when £ = 0.2, if the 
uncertainties take the form of = e(t)Ir (£ = x or y), we can improve the measurement 
period to T = = 1.049 using (fl"4l) . It is clear that T ^> T\ in these two cases. For some 
practical quantum systems such as spin systems in NMR, we can use strong control actions 
(e.g., K = 10 5 ) to drive the system from 1 1) into @ within a short time period T\ [8]. These facts 
make the assumption of no uncertainties in the control process reasonable. Moreover, the fact 
that the measurement period T is much greater than the control time required to go to |0) from 
|1) indicates the possibility of realizing such a periodic measurement on a practical quantum 
system. 

Remark 2: In the process of designing the control law for driving the system's state from 
|1) to |0), we employ an approach based on the Lyapunov methodology. An advantage of such 
an approach is that it is relatively easy to find a control law by simulation. It is worth noting 
that most existing applications of the Lyapunov methodology to quantum systems do not involve 
measurement. Here, we combine the Lyapunov-based control and projective measurements for 



D 01 02 03 04 05 06 07 

t 



Fig. 3. The control value u(t). 

controlling quantum systems, which in some applications make our method more useful than the 
Lyapunov-based control for quantum systems proposed in previous papers. In iHTTl . an approach 
based on time-optimal control design has also been proposed to complete this task. The advantage 
of such an approach is that we take the shortest time to complete the control task. However, 
it is generally difficult to find a complete time-optimal solution for high-dimensional quantum 
systems. For the above simple task, it has been proven that the time-optimal control employs 
a bang -bang control strategy |[T4|. Using the method in fl"4j| . we should take u = — 100 in 
t E [0, 0.016] and then use u = 100 in t E (0.016, 0.030]. In this case, the total time required is 
T[ = 0.030 « T x = 0.060). 

Remark 3: In the process of designing the Lyapunov control for driving the system's state 
from |1) to |0), we ignore possible uncertainties. By simulation, we find that small uncertainties 
can also be tolerated in this process. For example, if £ = 0.02 and the uncertainty e(t) is the 
noise with a uniform distribution on the interval [—0.02,0.02], the probability curves of |0) are 
shown in Fig. |4] when we apply the control obtained from Fig. [3] to the quantum system. The 
probabilities of |0) for the cases with uncertainties are very close to the probability of |0) for 
the case without uncertainties. By more simulation, we find that the final probability of |0) is 
(99.00 ±0.01)% for e{t)I x (\e{t)\ < e where £ = 0.02 or 0.2), the final probability of |0) is 
(99.00±0.02)% for e(t)I y (|e(f)| < 0.02) and the final probability of |0) is (99.00±0.13)% 
for e{t)I y (|e(?)| < 0.2). If we use a smaller probability of failure po (e.g., po = 0.5/?o) a s 




(a) Comparison between the case without noise and the case with y-axis noise 




without noise 
x-axis noise 
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(b) Comparison between the case without noise and the case with x-axis noise 

Fig. 4. The probability curves of |0) for the case without uncertainties (without noise) and the cases with uncertainties (noise) 
when we apply the control in Fig.[3]to the quantum system. The £-axis noise (S, =x or y) means the existence of e(t)h where 
e(f) is the noise with a uniform distribution on the interval [—0.02,0.02]. 



the terminal condition of the Lyapunov control or employ a bigger K for the same T\, these 
simulations suggest that it is possible to ensure that the Lyapunov control will drive the system's 
state into the sliding mode domain even when there exist small uncertainties. 

IV. Proof of Theorems 

This section will present the detailed proofs of Theorem 1 and Theorem 2. The proof of 
Theorem 1 involves the following steps: (I) Compare the probabilities of failure for H = I z + 
£ocos7o/ x + £osin7o/y and H = £ocos7o/ x + £osin7o/y (£o and Yo are constant); (II) Compare the 
probabilities of failure for H = el x and H = I z + e(t)I x (\e(t)\ < £); (III) Use the previous results to 
compare the probabilities of failure for H = el x and H = I Z + B x (t)I x + £ y (t)I y (yj E%(t) + Ey(t) < 
£); (IV) Use H = el x to estimate the measurement period. The basic steps for the proof of 
Theorem 2 include: (I) Formulate the problem of finding the "worst" case as an optimal control 



problem of minzf ; (II) Obtain the optimal control solution for nonsingular cases; (III) Exclude 
the possibility of singular cases; (IV) Use the "worst" case to estimate the measurement period. 
Considering that the arguments in the proof of Theorem 2 are useful for the proof of Theorem 
1, we will first present the proof of Theorem 2 and then prove Theorem 1. 



A. Proof of Theorem 2 

Proof: For H A = I z + e(t)I x , using p = -i[H A ,p] and ©, we have 

it xt-ijt \ _ I e ( f )yt -y t -ixt + ie(t)zt 

^x t + iy t -i t J y -y t + ix t -ie(t)zt -e(t)yt 

That is, 



(15) 



fi\ ( 

x t 



yit ) 



-1 

1 o 

e(f) 





-8(t) 





\ ( x \ 

Xf 



J 



yt 



(16) 



where (xo,y ,zo) = (0,0,1). 

We now consider e(t) as a control input and select the performance measure as 



(17) 



J(£)=Z f . 

From (fTTj) . we know that the "worst" case (i.e., the case maximizing the probability of fail 
ure) corresponds to minimizing Zf. Also, we introduce the Lagrange multiplier vector X{t) = 
(Ai (f ) , ^(0 ,^3(0 ) T and obtain the corresponding Hamiltonian function as follows: 

/ 

m(r(t),e(t),X(t),t)=X T (t) 



0-1 
1 -e(0 
e(t) 



\ ( x \ 

Xf 
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yt 



(18) 



where r(t) = (x t ,y t ,Zt). That is 

U(r(t),£(t)A(t)j) = -X l (t)y t + X 2 (t)x t + £(t)(h(t)y t -Mt)z,). (19) 

According to Pontryagin's minimum principle lloTI . a necessary condition for B*(t) to minimize 
7(e) is 

H(r*(f),e*(f),r(0,f) < H(r*(0,e(0,r(0,0- (20) 

The necessary condition provides a relationship to determine the optimal control £*(t). If there 
exists a time interval [t\,t2] of finite duration during which the necessary condition d20l) provides 



no information about the relationship between r*(t), £*(t), A*(f), we call the interval [?i,?2] 
a singular interval f[6TTl . If we do not consider singular cases (i.e., X^{t)y t — X%ii)zt = 0), the 
optimal control £*(t) should be chosen as follows: 



e*{t) = -£sgn(X 3 {t)y t -X 2 {t)z t ). 



(21) 



That is, the optimal control strategy for e(t) is bang-bang control; i.e., £*(t) = £ = +£ or — e. 
Now we consider H B = I z + £l x which leads to the state equation 
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yt 



(22) 



where (xo,yo,zo) = (0,0, 1). The corresponding solution is 
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yt 



1+e 
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Vl+e : 



sin cot 



\ T- 



■ cos cot + 



(23) 



l+e 2 / 



where CO = y/l +£ 2 . From (|23T) . we know that it is a monotonically decreasing function in t 
when t G [0, e2 \ - Hence, we only consider the case t G [0,?/] where tf G [0, £2 \ - 

Now consider the optimal control problem with a fixed final time tf and a free final state 
r f = ( x fiyf' z f)- According to Pontryagin's minimum principle, X*(tf) = J^r*(fy). From this, it 
is straightforward to verify that ^litf), Ag(?/)) = (0,0,1). Now let us consider another 

necessary condition X(t) = — Be1£M^)iM1i^ which leads to the following relationships: 



' hit) ^ 
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/n 1 o \ / 
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Ai(0 N 
A 2 (0 

^ h{t) j 



hit 

\ hit) j 

where (Ai(?y), A2(?/), hit/)) = (0,0,1). The corresponding solution is 

A 2 (r 
A 3 (r 



(24) 



V 



(25) 



We obtain 

h(t)yt - h,(t)z t = [sin (0t + £ 2 sin ©?/ + sin ©(/ /•-?)]. (26) 

It is easy to show that the quantity (A3 (/"))'/ — A 2 (?)zf) occurring in (I2TT) does not change its sign 
when tf £ [0,^=5] and £ £ [0,f/]. 

Now we further exclude the possibility that there exists a singular case. Suppose that there 
exists a singular interval [to,t\] (where to > 0) such that when t £ [?o,?i] 

^)=A 3 (0^-A 2 (0^ = 0. (27) 
We also have the following relationship 

h(t) = h(t)x t -?n(t)zt = (28) 
where we have used ( [TBI and the following costate equation 
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If to = 0, we have (xo,yo,Zo) = (0,0, 1). By the principle of optimality ED, we may consider 
the case t f = t\. Using (1271 (1281) and (Ai(fi), A 2 (?i), A 3 (?i)) = (0,0,1), we have x h = and 
y t} = 0. Using the relationship of xj +yj + zj = 1, we obtain z ?1 = 1 or — 1. If zt x = 1, the initial 
state and the final state are the same state |0). However, if we use the control e(t) = £, from 
(1231) we have zt x (e) = j^p cos ^1 + j^i <Zt x = l- Hence, this contradicts the fact that we are 
considering the optimal case minz/. If z fl = — 1, there exists < t\ < t\ such that zz = 0. By 
the principle of optimality lloHl . we may consider the case tf = t\. From the two equations (ITTT) 
and (1281) . we know that z- = 1 which contradicts Zf. = 0. Hence, no singular condition can exist 

fj 1 

if t = 0. 

If to > 0, using (T2TI) we must select e(t) = £ when ? £ [0,?q]. From (|26|) . we know that there 
exist no ?o £ (0,*/) satisfying A3 (?()).)% — A 2 (?o)z/ = 0. Hence, there exist no singular cases for 
our problem. 

From the above analysis, e(t) = £ is the optimal control when t £ [0, -y=^ \. Hence zf = 
Zf(e(0) > Zt(e) = zf. From (fTTT) . it is clear that the probabilities of failure satisfy pf = < 
pf = That is, the probability of failure pf is not greater than pf for t £ [0, , z J . When 



t e [0,-7=5], zf is monotonically decreasing and pf is monotonically increasing. When t = 
, n a , using ((231 we have zf = Irfj. That is, the probability of failure // = tttj- Hence, we 
can design the measurement period T using the case of H when < po < j^?z- 
Using £□]) and ((23]), for f G [0, ^=5] we obtain the probability of failure 

□ e 2 1— cos cot 

p " = \ 2 ' (30) 

Hence, we can design the maximum measurement period as follows 

^ (2) _ arccos[l-2(l + ^)p ] 
Vl+e 2 

For //a = £ (0^y (where \s(t)\ < £), we can obtain the same conclusion as that in the case 
H A = e(t)I x (where \e(t) \ < £). ■ 



B. Proof of Theorem 1 

To prove Theorem 1, we first prove two lemmas (Lemma [3] and Lemma [5]). 

Lemma 3: For a two-level quantum system with the initial state (xo,yo,zo) = (0,0, 1) (i.e., 
|0)), the system evolves to (xf,yf,zf) and (xf,yf,zf) under the action of H A = I z + £0 cos yiyl x + 
£0 sin yo/y (£0 is a nonzero constant) and H B = £0 cos yo/x + Eo sin yo/y, respectively. For arbitrary 

te[Q^,zf>zf. 

Proof: For the system with Hamiltonian H A =I z + £q cos JqI x + £0 sin yo/y, using p = —i [//, p] 
and ©, we obtain the following state equations 

/ 



yf 
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— 1 £o sin yo 

1 — eocosyo 
-e siny e cosy 

where (xq^^Zq) = (0,0, 1). The corresponding solution is as follows 

/ epsinyo 



yf 



(32) 



where coq 
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-«sm^-^cosa V + ^ 
V 1+£ o + 1+e o 



(33) 



For the system with Hamiltonian H B = 8ocos7o4 + eosin)t)/y, using p = —i[H, p] and ©, we 
obtain the following state equations 



£o sin 70 

-£ cos7o 

I -Eosinyo £ cos7o J 



4 

yf 



where (xq^q^Zq) = (OjOj !)• We can obtain the corresponding solution as 



37 



/ 
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/ 



(34) 



sin 7o sin £o? 

-cos7osin£of • (35) 
cos e t 

Since cos£o? = cos(— £ot), we may first assume £q > 0. We define F(t) and f(t) as follows. 

F(t)=z*-zf, (36) 



/(f) = F'(f) = Gosin£bf- 



Now we consider ? G [0, — t=^=j], and obtain 

V 1+e o 

/(f) = £o(cos£ f-cosc0o0 = 2£q sin 



l+£ 2 



:Sinft)of. 



(37) 



coo + £0 ^ - £0 
-f sin- 



f >0 



(38) 



2 2 

It is clear that /'(/) = only when t = 0. Hence f(t) is a monotonically increasing function and 



min/(0=/(0) = 0. 



Hence, we have 



(39) 



fit) > 0. 

From this, it is clear that F(t) is a monotonically increasing function and 

minF(f) =F(0) = 0. 

Hence F(f) > when t G [0, 7=^=7 ]- Moreover, it is clear that min f z^ = 7^|>_ » and 

V 1+e o 



V 1+£ o 

zf = cos£of is a monotonically decreasing function when t G [0, tJM. It is easy to obtain the 



following relationship 



k l*e[- 



2=«\ <Z?\ t= *^ <4\t-- 



'l+e.i 



Hence we can conclude that zf > zf for arbitrary t G [0, y|^] 



Let 70 = 0, and we have the following corollary. 

Corollary 4: For a two-level quantum system with the initial state (xo,yo,zo) = (0,0,1) 
(i.e., |0)), the system evolves to (xf,yf,zf) and (xf,yf,zf) under the action of H A = I z + £qI x 
(where £o is a nonzero constant) and H B = £qI x , respectively. For arbitrary t G [0, ^M, zf > zf ■ 
We now present another lemma. 

Lemma 5: For a two-level quantum system with the initial state (xo,yo,zo) = (0,0, 1) (i.e., 
|0)), the system evolves to (xf,yf,zf) and (xf,yf,zf) under the action of H A =I z + e(t)I x (where 
|e(f) | < e) and H B = el x , respectively. For arbitrary t G [0, §], zf > zf. 

Proof: First, we take an arbitrary evolution state (except |1)) starting from |0) as a new 
initial state. For H B = el x , the initial state can be represented as {x^y'^zo) = (0, — sin0o,cos 9q), 
where 6q G [0, 7t). We have 

/ 



(x*\ 

yf 



\ 









\ 


M 










-£ 


yf 


(40) 





e 


) 


\zf) 





The corresponding solution is as follows: 

(xf\ ( 
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(41) 



For H A = L + e(t)I x , we have 
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(42) 



where (x^,)^^) = (xo,yo,zo) = (sin0ocos<po,sin0 o sin<po,cos0o) and (po G [0,27!;]. From (142b . 
we have 



Zt\t=Q 



Z^ 

lim At . Z ° = £(0)y = e(0) sin0 o sin<p o . 



Af^O At 

From (143b and (|4T|) . it is easy to obtain the following relationship: 

4t ~ z At =zo + e (0) sin O sin cpoAt -z [l- ^^-} + sin G eAt + 0((At) 2 ) 
= At sin (£ + e(0) sin (po) + 0((At) 2 ) . 



(43) 



(44) 



When 6q e [0,7t), sin Gq > 0. Moreover, it is always true that £ + £ (0) sin <p > 0. If Go ^ and 
£ + £(0) sin <po 7^ 0, we have 

4>4r (45) 

If 0o = 0, we have (V ,y ,zo) = (xo,yo,zo) = (0,0, 1). According to Corollary @] and the proof 
of Theorem 2, we have 

4t>dt- (46) 

For £ + £(0) sin<po = 0, it corresponds to two cases: (a) £(0) = £ and <po = (b) £(0) = — £ 
and (pQ = f . In the following, we consider the case (a) (for the case (b) we have the same 
conclusion as the case (a)). Since <po = 3 j-, (xo,yo,zo) = (0, — sin0o,cos $o)- Using a similar 
argument to the proof of Theorem 2, we know that for H A = I z + e(t)I z with (xo,yo,zo) = 
(0, — sin0o ; cos0o)> tne °pti ma l control for the performance index J(e) = z/ takes a form of 
bang-bang control £(/) = £ = £ or — £. So we only need to consider a bang-bang strategy. 

For such a bang -bang strategy as H A = I z + el x , we have 

( % \ ( o -l o \ f 4^ 

yf = 1 -£ yf 
K zf ) [o e ) [z? / 



(47) 
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1 -i 
£ 

where (x^,y$,z$) = (xo,yo,zo) = (0, — sin O , cos Go) and <po e [0,2^]. We can obtain the corre 
sponding solution as 



yf 



( ^f^cos^ + ^^sin^+^a \ 



— ecos 61 

Vi+f 



" sin to? — sin 0q cos ft)? 



e-cosft) £ sin 0n •„ ,w 1 cosfti 



(48) 



Now, we consider the limit as At — > and obtain 

A _S _ £ 2 cos fo r, (l+£ 2 )(Af) 2 1 £sine / 1 j n 2\ t 
z At z At — i +£ i L 1 2 J Vl+£ 2 

+ £2^ _ cos e Q [i _ £W] + sin QeAt _ ^(a?)3 sin G + 0((At) 4 ) 
= AtsmG (£-£)+smG {At) 3 (£ + ££ 2 -£ 3 )+O((At) 4 ) 
>0. 

Hence, for arbitrary zo = cosGq (Gq G [0, 7t)), we have 



(49) 



(50) 



For t = ~, zf = — 1. Hence the relationship zf > zf is always true. 

We now define g(t) = — zf and assume that there exist t = t\ £ [0, ~) such that z^ < z^. That 
is, g(fi) < 0. Since g(t) is continuous in ? and g(0) = 0, there exists a time f* = sup{?|0 < t < 
h,g(t) = 0} satisfying g(t) < for t £ (?* , ?i] . However, we have proven that for any zf = zf and 
A? — > 0, zf +Af > zf +At , which contradicts g(Y) < for t £ (f*,*i]. Hence, we have the following 
relationship for t £ [0, |] 

2?>zf. (51) 



(52) 



Now we can prove Theorem 1 using Lemma [3] and Lemma [51 

Proof: For H A = I z + e x (t)I x + £ y (t)I y , using p = —i[H,p] and ©, we can obtain the 
following state equations 

/ & \ ( o -l e y (t) \ ( ^ 

1 -e x (f) 
v -ey(f) e x (0 y 

where (4^,4) = ( ,°: !)■ 

Define e(f) = Je$(t) + £j*(f) and £ x (r) = e(?)cosy, , 6y(r) = e(?)sin#. This leads to the 
following equation 

-1 e(f)siny* \ ( xf » 

1 -£(f)cosy 
, — £(j)siny* e(f)cosy f 



/ 



(53) 



where (x$,y^,z$) = (sinOocos <po, sin$o sin<po,cos 0q) and <po £ [0,2tt]. From (l53l) . we have 



Z;|;=0 = lim 



£ (0) cos y yo - £ (0) sin y x 



A? 

From (l54l) and (|4TI) . it is easy to obtain the following relationship: 



(54) 



4f - 4 = zo + e (0) cos y y A? - e (0) sin y jc A? 



-zo[l 



£ 2 (Af) 2 



-j + sin0 o £Af + O((A?) 2 ) 
Ar sin O (e + e(0) sin(<po - yo)) + 0((Ar) 2 )- 



(55) 



When O £ [0,tt), sin0o > 0- Moreover, it is always true that e + e(0) sin(<po — yo) > 0. If 
O 7^ and £ + e(0) sm((p - Yo) ^ 0, we have 

4>4- (56) 



If Oo = 0, we have (x^^y^zo) = (xo,yo,zo) = (0,0, 1). According to Lemma[3]and Pontryagin's 
minimum principle, we have 

4>4- (57) 

For £ + £(0)sin(<p - To) = 0, it must be true that e(0) = £ or -£. We consider £(0) = £ 
(for £(0) = — £ we have the same conclusion as £(0) = £). Moreover, we have Jo = <Po + § or 
7o = <Po-f • 

For 70 = (po + f , we first employ the fact that z t is independent on (po for H = e s'm(poI x — 
ecoscpoly and (xo,yo,zo) = (sin0ocos$o,sin0osin^,cos0o) since z t = cos(0o — £f). Then, it is 
easy to prove the relationship zfL > z^ using a similar argument to that in the proof of Lemma 
[3j For yo = <po — we have the same conclusion as in the case 7o = <Po + f • Thus, we obtain 

4t>4f (58) 

Now using a similar argument to that in the proof of Lemma [5l for arbitrary t G [0,| ), we 
have 

4 > zf- (59) 

From (fTTT) . it is clear that the probabilities of failure satisfy pf = < pf = ~x L - That 
is, the probability of failure pf is not greater than pf for t G [0, |). Hence, we can design the 
measurement period T using the case of H B . 

Using (fTTT) and (|4T1) . for ? G [0, -|), we obtain the probability of failure 

d 1— coset 

pf = j • (60) 

Hence, we can design the maximum measurement period as follows 

r (i) _ arccos(l-2p ) 



V. CONCLUSIONS 

This paper proposes a variable structure control scheme with sliding modes for the robust 
control of two-level quantum systems where an eigenstate is identified as a sliding mode. We 
present a design method for the control laws based on a Lyapunov methodology and periodic 
projective measurements to drive and maintain this system's state in the sliding mode domain. 



The key task of the control problem is converted into a problem of designing the Lyapunov 
functions and the measurement period. The Lyapunov function can be constructed to define a 
control law. By using simulation, we obtain an open-loop control to drive the controlled quantum 
system's state into the sliding mode domain. For different situations of the uncertainties in the 
system Hamiltonian, we give two approaches to design the measurement period, which guarantees 
control performance in the presence of the uncertainties. This sliding mode control scheme 
provides a robust quantum engineering strategy for controller design for quantum systems and 
has potential applications in state preparation, decoherence control, quantum error correction 
pT|. etc. Future work which can be carried out in this area is listed as follows. 1) The physical 
implementation of the proposed method on specific quantum systems. For example, spin systems 
in NMR (see, e.g., JH1, Il62l0 may be a suitable candidate to test the proposed approach. 2) The 
extension from two-level systems to multi-level quantum systems: The basic idea of sliding mode 
control can be extended in a straightforward way to multi-level quantum systems. However, it is 
much more difficult to obtain an analytical solution for a multi-level system. In [41], a numerical 
result has been obtained for a three-level quantum system to determine the measurement period 
and more complex systems are worth exploring by numerical methods. 3) The extension to 
dissipative quantum systems governed by the Lindblad equation [13] or described by a stochastic 
differential equation: For such cases, it is necessary to develop new methods to drive the system 
into the sliding mode domain since the Lyapunov-based control approach does not usually work 
|[63l . 4) The exploration of practical applications for the proposed approaches: The sliding mode 
may correspond to an eigenstate or a state subspace and the sliding mode design approach 
could be used in quantum state preparation and protection of encoded quantum information in 
a subspace. 



Appendix: Proof of > 



Proof: Take po as the variable and define 



F(p )=T^-T^ 



(62) 



For po E (0, jrgz), we have the following relationship 

f(po) = F'( P0 ) = - 1 = 
Je 2 p -(l+£ 2 )p 2 




1 



(63) 



It is clear from ([62]) and ([63]) that /(p ) > for p G (0, y^r) and F(p -> 0+) = 0. Hence 
F(p )>0foT Poe(0,j^). 

2 

When p = p' = j^, 

arccosfj-^v) 

r«(//) = (64) 

r (2 V) = ^=^- (65) 
V 1 + e 2 



Let x = -j—^ and 



l+e 
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G(£) = -JL=-arccos(|— ^). (66) 



We have 



K 



G(x) = —^y/l — x — arccosx. (67) 
V2 

For x G [— 1, 1], G(jc) is continuous in x and we also know that G(x) = only when x = ±1. 
It is easy to check G(x = 0) > 0. Hence, we know that for x G [— 1, 1], G(jc) > 0. That is, for 
e > 0, G(e) > 0. From the relationship G(e) > 0, we know T^(p') > T^(p') for £ > 0. 
Hence we concluded that for arbitrary po G (0, •jt^z]> 7^ > ■ 
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